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We present Spitzer IRAC (3.6, 4.5, 5.8, 8.0 /xm) and MIPS (24;um) obser- 

G> \ vations of the SCUBA submillinieter sources and /uJy VLA radio sources in a 

^1 5'x5' area in the Lockman Hole East region. Out of the ~40 SCUBA/VLA 

^ \ sources in the field, Spitzer counterparts were detected for nearly all except for 

t:;^- \ the few low-weight SCUBA detections. We show that the majority (80-90%) 

^ \ of the detected sources are cold (i.e., starburst-like) infrared- luminous galaxies 

^5 ' {LiR > 10^^ L0) at redshift 0.5 < 2; < 3.5, whose star-formation rate density 

' I (SFRD) is comparable to that of the optically-selected star-forming galaxies. 
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1. INTRODUCTION 

A major scientific goal of tlie Spitzer Space Telescope is to probe tlie evolution of infrared- 
luminous galaxies {Lm > 10^^ L0) at 2; > 1. Neither IRAS nor ISO provided enough 
sensitivity to detect a significant number of galaxies at 2; > 1, except for a small number 
of extreme objects. The much improved sensitivity and spatial resolution of Spitzer in the 
mid/far-infrared have improved the situation dramatically. With this goal in mind, we have 
imaged a 5'x5' area in the Lockman Hole East region with InfraRed Array Camera (IRAC; 
Fazio et al. 2004) and Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004). 
The Lockman Hole was chosen because of the low infrared cirrus background as well as the 
wealth of ancillary data. 

Up to now, z > 1 infrared luminous galaxies have been probed mainly through submil- 
limeter (e.g.. Small, Ivison, & Blain 1997) or radio (e.g., Barger, Cowie, & Richards 2000) 
observations. The former directly detects the emission by dust that is responsible for the 
infrared luminosity, while the latter detects the synchrotron emission (presumably from su- 
pernova remnants), which is known to correlate with the infrared luminosity (e.g., Helou, 
Soifer, & Rowan- Robinson 1985). These two types of selection provide complementary views 
in that submillimeter selection is biased toward ultra-luminous {Ljr > 3— 4x lO^^L©) sources 
at high redshift (with a median redshift of 2.4 (Chapman et al. 2003)), while radio selection 
is biased toward less luminous sources at intermediate redshift {z = 1 — 2). 

In this paper, we conduct the first analysis of SCUBA submillimeter and VLA radio 
sources based on the new Spitzer data. Companion papers also discuss stacking analysis 
of faint SCUBA sources (Serjeant et al. 2004), MAMBO millimeter sources (Ivison et al. 
2004a), IRAC sources (Huang et al. 2004), X-ray sources (Alonso-Herrero et al. 2004), 24 
yum-selected sources (Le Floc'h et al. 2004), and extremely red objects (Wilson et al. 2004) 
in the same field. 

Throughout the paper, Hq = 70 km s~^ Mpc~^, Qm = 0.3, and ^a = 0.7 were assumed. 



2. Data and Sample 

2.1. Observations and Data Reduction 

The details of the IRAC observations are presented by Huang et al. (2004), so only the 
MIPS observations are described here. The MIPS 24 //m images of the Lockman Hole East 
region were taken on UT 2003 November 30 (Program ID: 1077; AOR ID: 6619904). The 
large-source photometry mode was used to minimize the dither step lengths and maximize the 
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areal coverage overlap within the 5'4 x 5'A MIPS 24 /im field of view. The total integration 
time was 300 seconds per pixel. The 24 /xm imaging data were reduced and combined with 
the Data Analysis Tool (DAT) developed by the MIPS instrument team (Gordon et al. 2004). 
The final reduced 24 yum image is shown in Figure la while the three-color image produced 
from the IRAC 3.6, 8, and MIPS 24 fim. is shown in Figure lb. The 24 ^m sources were 
extracted from the image as described in Papovich et al. (2004). The detection limit in the 
24 ^m image is ~ 120 /iJy (3 0"). A preliminary IRAC/MIPS source catalog of this field will 
be released at http://mips.as.arizona.edu/cs/ero. 

A variety of other imaging data were also assembled (Le Floc'h et al. 2004). This ancil- 
lary data set includes the F606W HST/WFPC2 images as well as the UBVRIJHKs images 
obtained with a variety of ground-based telescopes. It also includes a deep (5 /xJy/beam) 
high-resolution (1'.'4 FWHM) VLA map (Ivison et al. 2002) and a reprocessed SCUBA map 
(Mortier et al. 2004). 



2.2. The SCUBA/VLA Sources in the Field 

The field imaged by both IRAC and MIPS contains 10 SCUBA sources in the 850 
Hm catalog of Scott et al. (2002), LE850.1/4/7/8/10/14/18/23/24/35i (Figure 1). The 
first seven sources are considered secure while the last three are regarded as marginal, the 
distinction being whether the signal-to-noise ratio is above or below 3.5. Among the seven 
secure sources, Ivison et al. (2002) detected radio counterparts for five (LE850. 1/7/8/14/18) 
and possibly for one more (LE850.4) with a lower significance. Our re-examination of the 
VLA map shows that LE850.35 also has a VLA source within a radius of 8", the region of 
95 % positional confidence according to Ivison et al. (2002). We tentatively identify this 
VLA source as the counterpart for LE850.35 although we could not confirm the reality of 
this submillimeter source unambiguously even with the reprocessed SCUBA map. LE850.35 
makes the total of radio-detected SCUBA sources seven out of ten. The field also contains 
a total of 38 VLA sources above 20 yuJy (4cr peak fiux density) using the map of Ivison et 
al. (2002). 



^The source names are based on Ivison et al. (2002). The full names according to SIMBAD are of the 
form [SFD2002]LHE N. 
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3. RESULTS and DISCUSSION 

3.1. Identification of SCUBA Sources 

Figure 2 shows the Spitzeras well as VLA and archived Subaru/Suprime-Cam (Miyazaki 
et al. 2002) i?-band images of all ten SCUBA sources in the field. The seven sources with 
radio detections have all been detected by IRAC while six (excepting only LE850.4) have 
been detected by MIPS. The IRAC/MIPS fiux densities of these sources are listed in Table 1 
together with the previously published 850 /zm and 20cm fiux densities. The following three 
sources have more than one IRAC source within a radius of 8": 

LE850.7 — The IRAC 3.6 and 4.5 yum images show extended emission north-east from 
the radio position, which is seen to be a nearby source in the i?-band image. However, the 
longer-wavelength images do not show emission in the same direction, indicating that this 
neighboring source does not contribute significantly above ~ 5 ^m. 

LE850.8 — Three IRAC sources (8a/8b/8c) are seen within a radius of 8". The 24 
ixvn signal is clearly extended in a way consistent with all three sources contributing. The 
strongest radio source corresponds to component 8b. There is also faint (3 a) radio emission 
at the position of component 8a, which is the brightest source at 24 ^m. Component 8a was 
thought to be the counterpart for the ROSAT X-ray source with a spectroscopic redshift of 
0.974 (Lehmann et al. 2001), but the later XMM-Newton observation showed that the X-ray 
source is actually coincident with component 8b (Ivison et al. 2002). This leaves a possibility 
that some of the radio fiux from component 8b may be of AGN origin. The /?-band image 
shows several sources at the position of component 8b. 

LE850.14 — The two radio counterparts (14a/ 14b) have both been detected by IRAC. 
The 24 /im emission is centered between the two IRAC sources. Considering the equally 
strong radio fiuxes, both components are expected to contribute significantly in the submil- 
limeter. 

In summary, we have 7 secure SCUBA sources with a total of 9 radio components 
(LE850.1/4/7/8a/8b/14a/14b/18/35), all of which have S])iteer counterparts. Spitzer coun- 
terparts have been found only for the sources with radio detections. The three sources which 
are not detected by VLA and Spitzer (LE850. 10/23/24) are all weak in the reprocessed 850 
/im map (Figure 2). 
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3.2. Spectral Energy Distribution and Redshifts 

For the following statistical analysis, we restrict our discussion to the VLA 20 fiJj 
sample, which contains 7 of the 9 SCUBA radio components (LE850.1/7/8b/14a/14b/18/35). 
There is one VLA source which is relatively bright in the radio (230 ;uJy) but undetected in 
most other bands including the IRAC 5.8 and 8.0 yum (this is the only source which has not 
been detected in all four IRAC bands). We exclude this source; The VLA sample consists of 
37 objects, of which all are detected in the four IRAC bands, and 29 are detected at 24 /im. 

The SEDs of the VLA sources can be classified into two types: those showing a clear 
near-infrared stellar continuum hump around 1.6 //m (32 sources, 86%) and those showing a 
featureless power-law continuum (5 sources, 14%). Photometric redshifts can be derived for 
the former using the near-infrared hump (e.g., Sawicki 2002; Le Floc'h et al. 2004); For the 
latter, the redshifts are mainly constrained by the submillimeter and radio points, but the 
estimates are significantly more uncertain. 

Figure 3 shows the composite SEDs of these two types. The SED of the sources with a 
near-infrared continuum hump is similar to that of a "cold" ultra-luminous infrared galaxy 
(ULIRG) like Arp 220 while the SED of the power-law continuum sources resembles that 
of a "warm" ULIRG like Mrk 231. This cold/warm classification often used with ULIRGs 
roughly separates starburst-dominated galaxies (cold) and AGN-dominated ones (warm) by 
the appearance of the infrared SEDs. Among the 32 cold-type sources, there are 8 sources 
(22% of the 37 VLA sample) which were not detected at 24 jum, but their restframe near- 
infrared SEDs are similar to those of the 24 yum detected sources. Considering that these 
sources have expected 24 /xm fluxes (based on the Arp 220 SED) close to or below our 
detection limit, we include them in the cold category. 

As Figure 3a shows, the cold-type sources show remarkably similar SEDs from the 
restframe visible to near- infrared that are well-fltted by the SED of Arp 220. Such unifor- 
mity is in contrast to the diversity of the SEDs seen with the X-ray and 24 //m selected 
sources (Alonso-Herrero et al. 2004; Le Floc'h et al. 2004). Among the SCUBA sources, 
LE850.1/7/14a/14b/35 have this type of SEDs. (LE850.4 is also of this type, but it is not 
formally in the radio-selected sample.) With the SCUBA sources, the good flt to the Arp 220 
SED extends to submillimeter and radio while most of the VLA-only sources have signifl- 
cantly lower radio fluxes, which reflects the fact that they are not as infrared-luminous as 
Arp 220. Figure 3a also shows that higher redshift sources (i.e., those with radio points at 
shorter wavelengths) have radio SED points closer to that of Arp 220. 

The small SED dispersion with the cold-type sources demonstrates not only the sim- 
ilarity of the SEDs but also the accuracy of the photometric redshifts based on the near- 



-6- 



infrared continuum hump. Table 1 shows that the redshifts derived this way for the SCUBA 
sources are consistent with those based on the submm/radio flux ratio (Aretxaga et al. 
2003). LE850.14a and 14b show very similar SEDs, which indicates that their redshifts are 
both z ~ 2.5. In the case of LE850.14a, this was confirmed by the spectroscopic redshift of 
2.38 (Ivison et al. 2004b). This suggests a possibility that these two sources are physically 
associated. 

If we regard the warm sources as AGN galaxies, the AGN fraction of the 20 /xJy radio 
sample is 14% (5 out of 37). Among the 9 radio components associated with the SCUBA 
sources, the AGN fraction is ~ 30% (3 out of 9) as shown in Table 1. A comparison with the 
150 ksec XMM data (Alonso-Herrero et al. 2004) shows that only 2 out of the 5 power-law 
continuum sources are detected, and that there are three cold sources detected in the X-ray. 

Figure 4a shows the redshift distribution of the 32 cold sources. Assuming that these 
are starburst-dominated, the figure can be regarded as the redshift distribution of infrared- 
luminous star-forming galaxies selected in the radio. The radio- infrared correlation suggests 
that the infrared luminosities range from 10^^ to 10^^ Lq, except for a few sources at ^ < 0.5, 
which are below 10^^ Lq. The SCUBA sample constitutes only a small subset of the radio 
sample at high redshift. This is because (1) due to the relatively high detection limit [Lir > 
3 — 4 X 10^^ Lq), the submillimeter observation misses lower luminosity sources which radio 
observation can detect at z < 3 (e.g., Barger, Cowie, & Richards 2000), and (2) because of 
the strong luminosity evolution of the infrared-luminous galaxies, large luminosity systems 
detectable above the submillimeter detection limit are more abundant at higher redshift 
(e.g.. Small, Ivison, & Blain 1997). 



3.3. Star Formation History 

Using the derived redshift distribution and converting the 1.4 GHz radio luminosity into 
a star formation rate with the conversion factor of Yun, Reddy, & Condon (2001), we have 
calculated the star formation rate density (SFRD) as a function of redshift (Figure 4b). A 
K-correction was applied based on a z/^°'^ spectrum. Our SFRD at 2; = 1 is very close to 
the value derived by Barger, Cowie, & Richards (2000) with a sample of //Jy VLA sources 
with spectroscopic redshifts. 

The observed SFRD of the VLA sample shows a small peak at ^ = 2, but otherwise stays 
fiat at 2; = 1 — 3 above the observed SFRD of optically selected star-forming galaxies. When 
corrected for dust extinction, the latter exceeds the former by a factor of a few. However, 
the SFRD of the infrared-luminous sample includes only the most extreme high-luminosity 
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objects at high redshift. A substantial completeness correction is still required, which has 
already been applied to the optical sample. 

To estimate the size of the completeness correction, the observed SFRD was compared 
with a model prediction by Lagache et al. (2004). This model incorporates a variety of the 
latest observational constraints including the MIPS 24 ^m number counts (Papovich et al. 
2004). However, the behavior of the model SFRD in Figure 4b is not unique to this model 
and is similar to what was derived previously (e.g., Franceschini et al. 2001). The observed 
points (solid circles) and the model-predicted SFRD adjusted for the 20 /xJy radio detection 
limit (solid line) are seen to agree atz = 2 — 3, especially when we include the 24 ^m detected 
infrared-luminous galaxies which should have been in the radio sample (based on the 24 fira 
^ IR ^ radio luminosity conversion) but are not (e.g., sources larger than a few arcseconds 
would have been resolved away in the radio map). This suggests that if the completeness 
correction is applied, the total SFRD of this infrared-luminous galaxy population, which is 
considered to be distinct from the optically selected one (e.g., Franceschini et al. 2001), is 
flat at 2; = 2 — 3 at a level comparable to the SFRD of the optically selected galaxies. At 
2; ~ 1, however, the model overpredicts the observed SFRD by almost a factor of two. We 
will examine in the future papers whether this discrepancy is particular to this fleld (e.g., 
cosmic variance) or persistent even with a larger sample. 



4. Summary 

Except for the few low-weight SCUBA detections, nearly all the SCUBA and VLA /iJy 
sources have been detected by Spitzer. The majority (86%) of the sources show remarkably 
similar SEDs well fltted by the SED of a cold ULIRG like Arp 220. This characteris- 
tic was exploited to derive the redshift distribution, which then enabled us to derive the 
star formation history of this galaxy population using the radio-infrared correlation. These 
submillimeter /radio selected galaxies are mostly star- forming (i.e., cold) infrared-luminous 
galaxies [Lj^i > 10^^ Lq) at 0.5 < z < 3.5, whose SFRD is probably comparable to that of 
the optically-selected star-forming galaxies. 

We would like to thank S. C. Chapman for communicating us the revised redshift of 
LE 850.18 before publication, and G. Neugebauer for commenting on the manuscript. This 
work is based in part on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA through Contract Number 
960785 issued by JPL/Caltech. 
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Fig. 1. — (a) MIPS 24 /iin image of the Lockman Hole East; (b) Three-color image of the 
same field produced from the IRAC 3.6 /im, 8.0 /xm, and MIPS 24 /zm images. The ten 
SCUBA sources (LE850. 1/4/7/8/10/14/18/23/24/35) are indicated with white circles with 
a radius of 8". The field of view is ~5'x5'. 
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Fig. 2. — Postage stamp images of the ten SCUBA sources. Shown here are (from left to 
right) VIA, IRAC 3.6, 4.5, 5.8, 8.0 ^m, MIPS 24 yum, and R-band images. The VIA images 
are from the data of Ivison et al. (2002) while the R-band images are from the data taken 
from the Subaru telescope archive. The 850 /im contour maps (Mortier et al. 2004) are 
overlaid on the VLA images. The contours start at 2.5(7 above the sky and increases by a 
factor of 1.2. The size of each image is 20"x20". The solid circle indicates the region of 95 
% positional confidence (r = 8") by Ivison et al. (2002). The small white circles indicate the 
positions of the radio sources seen in the VLA map. In all the images, north is up and east 
is left. 
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Fig. 3. — The SEDs of 37 VLA sources detected above 20 fiJj, which contain 7 SCUBA 
sources: (a) The composite SED of 32 sources with a near-infrared continuum hump. The 
SED is well fitted by that of a cold ULIRG Arp 220 (solid line; the model by Silva et al. 
(1998)); (b) The composite SED of 5 power-law continuum sources. The SED resembles 
that of a warm ULIRG Mrk 231 (solid line; the observed SED). The solid circles denote 
SCUBA sources while the open circles denote VLA sources without SCUBA detections. The 
SEDs are shifted to the restframe using photometric/spectroscopic redshifts, and normalized 
around 1.6 /xm (i.e., the peak of the stellar continuum) with the cold sources and at 1-10 /im 
with the warm sources. For the three warm VLA sources without submillimeter detections, 
for which redshift cannot be estimated, zero redshift was used. 
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Fig. 4. — (a) The redshift distribution of 32 cold-type VLA sources. The dark shaded region 
indicates the 5 SCUBA sources, (b) Star formation rate density (SFRD) of the cold-type 
galaxies plotted as a function of redshift with three redshift bins (0.5-1.5, 1.5-2.5, 2.5-3.5). 
The bins contain 17, 9, and 3 sources, respectively, and the resultant SFRDs are plotted 
as solid circles with the error bars based on the Poisson counting uncertainty (i.e., SFRD 
divided by vN). The open squares are from Barger, Cowie, & Richards (2000), and the open 
triangle is from Yun, Reddy, & Condon (2001). No completeness correction was made for the 
z > points. The lower shaded band indicates the range of observed SFRD with optically 
selected galaxies while the upper shaded band indicates the range with dust extinction 
correction (e.g., Steidel et al. 1999; Giavalisco et al. 2004; Bouwens et al. 2004). The dashed 
line indicates the total SFRD predicted by the model of Lagache et al. (2004) while the solid 
line shows the model SFRD measurable above the radio flux limit of 20 yuJy. The arrows 
indicate the increase of SFRD when we include the 24 /im selected infrared-luminous sources 
from Le Floc'h et al. (2004) which have been missed by the radio observation. 



Table 1. SPITZER flux measurements of the SCUBA sources 



Source 


3.6 ^m 


4.5 ^m 


5.8 ^m 


8.0 ^m 


24 yum 


850 Mm 


20 cm 




z •" 


Zphot 


Type 


(LE850) 


(MJy) 


(MJy) 


(MJy) 


(A*Jy) 


(MJy) 


(mJy) 


(A^Jy) 


(submm/radio) 




(this work) 




1 


3.8±0.6 


8.2±1.1 


11.8±2.5 


10.2±1.7 


193±33 


10.5±1.6 


73±10 


2.6±g:t 




2.6 


cold 


4 


7.7±1.0 


10.2±1.4 


11.5±2.5 


9.6±1.6 




8.3±1.8 


19±8 


4.7±J:i 




2.6 


cold 


7 


65.9±6.9 


77.0±8.1 


60.1±7.8 


52.4±5.9 


312±39 


8.1±1.9 


135±13 


3.6±?i 




1.8 


cold 


8b 


17.1±2.0 


22.3±2.6 


37.2±5.3 


89.5±9.6 


282±59 


5.1±1.3 


58±12 


3.7±J:f 




-3 


warm 


8a 


72.0±7.5 


72.6±7.6 


79.0±9.7 


112.0±11.9 


534±117 




22±11 




0.974 




warm 


8c 


26.7±3.0 


20.7±2.4 


17.2±3.1 


17.8±2.4 


161±47 










0.9 


cold 


14b 


9.8±1.3 


14.7±1.8 


22.2±3.7 


18.5±2.5 


166'= 


6.3^= 


72±12 


2.4±J:1 




2.5 


cold 


14a 


15.2±1.9 


20.2±2.4 


27.6±4.3 


18.6±2.5 


83^= 


3.2^= 


36±12 




2.38 


2.5 


cold 


18 


6.1±0.9 


5.8±0.9 


19.9±3.4 


26.4±3.3 


125±33 


4.5±1.3 


47±10 


3.0±l:l 


2.69 


-3 


warm 


35 


23.5±2.7 


25.9±2.9 


26.3±4.1 


34.8±4.2 


155±31 


6.7±2.3 


56±10 






3.0 


cold 



hl^ 



Note. — The photometric uncertainties include the 10% calibration errors. The 850 /im fluxes and 20 cm fluxes are from Scott et al. (2002) 
and Ivisonet al. (2002). 

'^Photometric redshifts (le2 type) from Aretxaga et al. (2003) 

•"Spectroscopic redshifts for LE850.8a (Lehmann et al. 2001), LE850.14a (Ivison et al. 2004b), and LE850.18 (Chapman et al. 2004). 

^The measured total flux densities of LE850.14 at 24 /im (249±39 /iJy) and 850 /im (9.5±2.8 mJy) were distributed based on the 20cm flux 
density ratio. 



